The December 2009 gamma-ray flare of 3C 454.3: the 
multifrequency campaign 
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ABSTRACT 

During the month of December, 2009 the blazar 3C 454.3 became the brightest 
gamma-ray source in the sky, reaching a peak flux F ~ 2000 x 10~®ph cm~^ 
for E > 100 MeV. Starting in November, 2009 intensive multifrequency 
campaigns monitored the 3C 454 gamma-ray outburst. Here we report the results 
of a 2-month campaign involving AGILE, INTEGRAL, Swift /XRT, Swift /BAT, 
RossiXTE for the high-energy observations, and Swift /UYOT, KANATA, GRT, 
REM for the near-IR/optical/UV data. The GASP/WEBT provided radio and 
additional optical data. 

We detected a long-term active emission phase lasting ~1 month at all wave- 
lengths: in the gamma-ray band, peak emission was reached on December 2-3, 
2009. Remarkably, this gamma-ray super-flare was not accompanied by corre- 
spondingly intense emission in the optical/UV band that reached a level substan- 
tially lower than the previous observations in 2007-2008. 

The lack of strong simultaneous optical brightening during the super-flare and 
the determination of the broad-band spectral evolution severely constrain the 
theoretical modelling. We find that the pre- and post-flare broad-band behavior 
can be explained by a one-zone model involving SSC plus external Compton 
emission from an accretion disk and a broad-line region. However, the spectra 
of the Dec. 2-3, 2009 super-flare and of the secondary peak emission on Dec. 9, 
2009 cannot be satisfactorily modelled by a simple one-zone model. An additional 
particle component is most likely active during these states. 

Subject headings: galaxies: active - galaxies: jets - galaxies: individual (3C 454.3) - 
radiation mechanisms: non-thermal 
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Introduction 



The flat spectrum radio quasar 3C 454.3 (at a redsliift z = 0.859) turns out to be 
among tlie most active blazars emitting a broad-spectrum ranging from radio to gamma-ray 
energies. Blazars are a sub-class of active galactic nuclei, with the relativistic jet aligned to 
the line of sight. Their spectral energy distributions (SED) typically show a double humped 
shape, with the low en ergy peak lying between r adio and X-rays, and the high energy peak 



in the GeV-TeV band lUrry and Padovanil (119951 ). Detailed description of blazar leptonic 



emission models can be found in Maraschi et al. (1992); Marsher & Bloom (1992); Sikora et 



al. (1994). The obseryed spe c tra can also b e mod e 



Mucke and Protheroel fl200ll ): 



Mucke et al 



(120031 ): 



led in the framework of hadronic models 



Bottchej (120071 ). 



Starting in 2004-2005, 3C 454.3 showed a long period of optical activity, with variability 
timescale ranging from several months to less than one day. In M ay 2005, the source 



Villata et al. 



(12006r i. 



reached a peak magnitude i? ~ 12 showing strong 1-day variability 
A radio peak was detected 9 months aft e r the optical peak, with a flux of 22 Jy at 37 



GHz, and 20 Jy at 43 GHz 



Villata et al 



(120071 ). The source was then quiesc ent from the 



begin ning of 2006 until mid-2007 (with an R magnitude between 1 5 and 16, iRaiteri et al. 



20081). Starting in the second half of 2007 



in a high g amma-ray stat e by AGILE 



Fermi-LAT 



Atwood et al. 



mm 



Vercellone et al. 



Tavani et al. 



Abdo et al. 



(120081 ) ■ has been detected 



( 120091 ) ■ and, subsequently also by 



( 120091 ). Typically, the level of gamma-ray 



activity (with a flux of 300-600x10 ^ ph cm ^ s ^ for E > 100 MeV) has been observed to 
be correlated with the optical emission. Relatively large gamn i a-ray fluxes were detected 



during the AGILE observations of 2007 



Vercellone et al. 



(12008 



20091) 



Donnarumma et al. 



( I2OO9I ). During 2008 , the source that wa s bright at the beginning of the year, started to 



fade in optical b and 



Vercellone et al. 



Villata et al. 



( I2OO9I ). AGILE detected the fading in gamma-rays too 



(120101 ). Fermi reported an averaged spectrum above 200 MeV with a 
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photon inde x a 



Abdoetal 



mm 



2.3 and a spectral break at Ec ~ 2.4 GeV, obtained in August, 2008 



2. The multifrequency campaign 



The intensive monitoring of garried out by our group cov ered the period of 



extraordinary gamma-ray activity in November-December, 2009 iStriani et aLl ( l2009l ). The 
campaign involved AGILE for the gamma-ray band, Swift/BAT, RossiXTE/HEXTE and 
INTEGRAL/IBIS in the hard X-ray band, RossiXTE/PCA and Swift /XKT in X-rays, 
Swift /UVOT in the optical and UV bands, the KANATA observatory and GRT in the 
optical, REM in the near infrared and optical. AGILE observed the source every day in 
spinning mode, scanning about 70% of the whole sky every 6 minutes. I NTEGRAL point ed 



the sou rce in response to a Target of Opportunity observations (ToO) iVercellone et al 



( l2009bl ). and observed it from 2009 December 6 until 2009 December 12. The RossiXTE 
satellite observed 3C 454.3 on 2009 December 5 and then daily from December 8 until 
December 17, 2009 for typical integrations of ~3 ks. Swift started to observe on 2009 
November 27, in response to a ToO, and pointed at the source every day (UVOT 
performed most of the observations with the UV filters). The KANATA 1.5 m telescope 
performed a long source monitoring in the V-band, with a time-step of 1 day. The fully 
automated 14" GRT (Goddard Robotic Telescope) performed observations in the V and R 
bands, quasi- simultaneously with Swift, starting on 2009 November 30. REM started the 
monitoring on 2009 December 10 in response to a ToO, and observed t he source every day 



in th e VRIJHK filters. The GLAST-AGILE Support Program (GASP 



Villata et al. 



2008 



20091 ) performed an intensive monitoring-campaign of the source in 2009-2010. We used a 
sub-sample of their data: Optical observations reported in this paper were performed at: 
Lulin, New Mexico Skies, Roque de los Muchachos (KVA), and St. Petersburg. GASP 
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radio data were taken at Mauna Kea (SMA, 230 GHz), Noto (43 GHz), Metsahovi (37 
GHz) , and UMRAO (4.8, 8.0, and 14.5 GHz). 



3. Data analysis 

AGILE/GRID data were analyzed using the Build-19 software and the response matrix 
vlO calibrated in the energy range 100-3000 MeV. Well reconstructed gamma-ray events 
were selected using the FM3.119 filter. All the events collected during the passage in the 
South-Atlantic Anomaly were rejected. We filtered out the Earth-albedo, rejecting photons 
coming from a circular region of radius 85 deg and centered on the Earth. We rejected 
photons coming from outside the 35 degrees from the optical axis. Gamma-ray data were 
analyzed with integrations of 1 or 2 days, depending on the so urce flux. We used the 



1996 



for the 



standard AGILE Maximum-Likelihood procedure (ALIKE) (see iMattox et al. 
concept deflnition) for each data set. The integration over 5 weeks from 2009-11-18 to 
2009-12-23 UTC yields a photon index 1.88 ± 0.08 (all the errors reported in the paper are 
at la, except where stated). 



The INTEGRAL-IBIS 



Ubertini et al 



(I2OO3I ) data were processed using the OSA 



software version 8.0. light curves (from 20 to 200 keV) and spectra (from 18 to 200 keV) 
were extracted for each individual science window of revolutions 873 and 874. 

The Swift-BAT survey data were obtained applying the "BAT FOV" option. The data 
have been processed by batsurvey script available through HEASOFT software package 
with a snapshot (single pointing) interval. To estimate the background, ten background 
points around the source in a radius of 50' are selected. The source, ten background points 
and the bright hard X-ray sources (for cleaning purpose) are included in the input catalog of 
batsurvey. The BAT count rate in the 14-195 keV band has been converted into the energy 
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flux assuming a power-law photon index 1.7 as determined from the INTEGRAL-ISGRI 
data, see below). In order to match with the HEXTE range, the BAT hard X-ray flux has 
been rescaled in Fig. [T]to the 20-40 keV band. 



RossiXTE-FCA Ijahoda et all (119961 ) and HEXTE iRothschild et all (119981 ) data were 
analyzed following the same procedure described in Vercellone et al. (2010). The data 
analysis was restricted to the PCU2 in the 3-20 keV energy range for the PGA and to 
the Gluster B in the 18-50 keV range for the HEXTE. The net exposure times were 27.3 
ks for PGA and 7.3 ks for HEXTE. The background subtracted source spectra obtained 
with both the instruments were simultaneously filo with an absorbed power-la w model. 



with the photoelectric absorption fixed to 0.134x10 cm 



Villata et al. 



(120061). After the 



introduction of a 2% systematic error, the best-fit value for the photon index is 1.74 ± 0.01. 

The Swift-XBT data were processed using the most recent calibration files available. 
We utilized Swift Software version 3.5, FTOOLS version 6.8, and XSPEG version 12.5. In. 
We fitted the data with an absorbed power-law model. We obtained photon indices between 
1.51 ± 0.09 and 1.73 ± 0.11, and excess absorption between (0.09 ± 0.06) xlO^^ and (0.17 
± 0.03) X 10^^ cm^^ (all uncertainties on xrt spectral fit are at 90% level). 

Swift-UVOT data from each observation sequence were processed by the standard 
UVOT tool uvotsource using the same version of the Swift software as for the XRT 
analysis. An extraction region of radius 5 arcsec centered on the source and a suitable 
backg r ound region were used. Magnitudes are based on the UVOT photometric system 



Poolel (l2008[ ). 



^The resulting fit is not statistically acceptable (Xred~-'--^3/54 d.o.f.); however, the addi- 
tion of a 2% systematic error to the data, which is well within the expected uncertainties in 
the spectral calibration, is sufficient to make the fit fully acceptable (Xred=l-02/54 d.o.f.). 
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The optical 



Watanabe et al. 



photom etry of the Kanata Observatory data was performed using TRISPEC 



(|2005[ ). The observations were pipehne-reduced, including bias removal 



and fiat-field corrections. We derived the V-band magnitu de from differential photom etry 



Gonzales-Perez et al. 



200l|). 



with a nearby reference star, USNOB 1061-0614254 ( see 
Photometric stability of this star has been confirmed by our simultaneous photometry for 
another nearby star, USNOB 1061-0614207. 



All REM raw opti c al an d NIR frames, obtained with ROSS iTosti et al.l ( 120041 ) and 



REMIR iConconi et al 



(120041 ) respectively, as well as images from GRT were corrected 



for dark, bias, and fiat field following standard recipes. Instrumental magnitudes were 



obtained via aperture photometry, and absolute call 



of secondary standard stars in the field 



j ration has been performed by means 



Raiteri et al. 



(119981). 



Even if it is not reported in Figure 1, the Metsahovi radio data at 37 GHz show a high 
flux with an increasing trend from 2009 December 1 until 2010 January 14, and a mean flux 
of ~ 20 Jy duringthe flrst week of December, 2009. The mean 230 GHz flux is ^ 25 Jy. 



These radio fluxei 



are comparable to the peak flux measured in 2006 



Villata et al. 



(l2007h 



All UV/optical/NIR data present ed here were correcte d for the Galactic extinction 



toward 3C 454.3 assuming Ay = 0.349 



Schlegel et al. 



(119981 ). 



4. Results 



The multifrequency light curves of §re reported in Figure HI The exceptional gamma- 
ray flaring activity is produced during an extended period lasting several weeks. The 



^The long term radio and optical light curves of 4uring the 2009-2010 observing season 
will be presented in a forthcoming paper (Raiteri et al., in preparation). 




Fig. 1. — Multifrequency data collected during the campaign (Upper panel:) Optical V 
data from: KANATA (diamonds), REM (star), CRT (squares), UVOT (triangles). (Middle 
panel:) X-ray data from XRT (0.3-10 keV flux, black triangles), PGA (3-20 keV, red squares), 
HEXTE (20-40 keV, green squares), BAT (14-195 keV rate, rescaled to the 20-40 keV band, 
cyan triangles), ISGRI (20-200 keV, blue diamonds). (Bottom panel:) gamma-ray data from 
the AGILE-GRID (0.1-3 GeV, black diamonds). 
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broad-band coverage is excellent, and for the first time allows a detailed mult i- wavelength 
study of a gamma-ray blazar peaking above the Vela pulsar flux (for c omparison, see the 



multifrequency coverage for the gamma-ray super-fiare of PKS 1622-29, 



Mattox et al. 



19971). 



The optical data show variability timescale as short as one day or less. Due to the 
typical interval between optical observations (1 day), we did not show in Fig. [T]a detailed 
representation of the variability. The optical V data indicate a low-state of the source (V 
> 15 mag, flux < 3.9 mJy) until MJD 55160, after which the flux started to increase. The 
optical flux increased of about 50% in less then one day from MJD 55166.4 until MJD 
55167.4, reaching the value V = 13.7 mag (e.g., 12.8 mJy), followed by a fast flux decrease 
to V" = 14.3 mag (e.g., 7.4 mJy) at MJD 55169.0. Another optical peak was reached at 
MJD 55172.5 (V=13.7 mag), and then was followed by a minimum at MJD 55175.0, with 
V=14.4 mag (e.g., 6.7 mJy). 

In general, the X-ray flux followjf] the rising part of the optical emission in the interval 
MJD 55150-55169. Starting on MJD 55169, the X-ray flux started to fade, with the optical 
emission remaining in a relatively high state for 4-5 days. Our INTEGRAL hard X-ray 
data sample the fading phase of the high-activity period. 

We performed a Discrete Correlation Function (DCF) analysis for the optical-gamma, 
and soft x-gamma data set. No significant delay is found within 2 days (compatible with 
the average bin of gamma-ray light curve). 

We focused on the time-dependent spectral analysis of this exceptional activity of 3C 
454.3 , and used simultaneous broad-band data to obtain a detailed account of the source 
variability. We obtained the SED for four periods. The first period (interval-1) is for a 



^We note that no X-ray data were obtained in exact correspondence with the gamma-ray 
super-flare of Dec. 2-3, 2009. 
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logv [Hz] 



T ' 1 1 1 1 1 ' 1 ' 1 ' r 

■ flare II 




10 12 14 16 18 20 22 24 26 

logv [Hz] 

Fig. 2. — Multifrequency simultaneous spectra of fit four different epochs superimposed 
with spectral modelling. Top panel: pre-flare (red open circles) and super-flare (black filled 
squares). Bottom panel: secondary flare (blue filled squares) and post-flare spectra (green 
open circles). Models details are explained in the text. 
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5-day integration of GRID data, centered on 2009-11-27 09:36 UT (MJD 55162.4), and 
using the Swift observation at MJD 55162.9, and quasi- simultaneous GASP optical data 
(from KVA and New Mexico Skies) obtained within 16 hours from the Swift observation 
(pre-flare SED). 

The second period (interval-2) was obtained for the gamma-ray super-fiare episode 
integrating the GRID data for 1 day centered at 2009-12-02 16:48 UT (MJD 55167.7), and 
using of the simultaneous Swift and GRT observations at MJD 55167.0 (first fiare SED). 
The Kanata observatory measured a fiux 40% higher in V band 10 hours after the GRT 
observation. 

The third SED was obtained (interval-3) integrating GRID data for 2 days centered at 
2009-12-06 16:48 UT (MJD 55172.7), to match the local maximum of the optical light curve 
with V=13.7 mag, that is apparently coincident with the secondary gamma-ray maximum 
near MJD 55174. For this interval we used the INTEGRAL-ISGRI data collected between 
MJD 55171.7 and MJD 55174.2, Swift data at MJD 55173.9, RossiXTE data at 55173.4, 
and GASP optical data (from Lulin and GRT) obtained within 26 hours from the Swift 
data. 

The last SED (interval-4) was obtained integrating GRID data for 5.5 days centered 
at 2009-12-15 18:00 UT (MJD 55180.8), and making use of the Swift observation at MJD 
55179.1, the RossiXTE observation at MJD 55179.2, the GASP optical data (from St. 
Petersburg and Lulin) simultaneous with Swift within 12 hours, and the near-infrared 
observations from REM at MJD 55181.0 (post fiare SED). 

Radio data for the SEDs were taken by the GASP, using observations from Mauna Kea 
(SMA, 230 GHz), Noto (43 GHz), Metsahovi (37 GHz) UMRAO (4.8, 8.0, 14.5 GHz), 
simultaneous within a few days with the XRT observations. Due to the slow variability of 
the radio data, we use interpolated radio values in the SEDs. The results are shown in Fig. |2j 
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5. Discussion 



The multifrequency data of the extensive campaign on §how a remarkable behavior 
of the source. Starting from MJD 55150, first an overall rise of the gamma-ray emission 
and then of the X-ray and optical fiuxes is detected. This rise culminates with peak 
optical/X-ray/gamma-ray emission during a 10-day period centered around MJD 55173. 
Subsequently, the overall fiux decreased and reached again a relative-minimum state around 
MJD 55200. During the 2-month period the optical and X-ray fiuxes vary within a factor of 
3, whereas the 7-ray fiux grows by a factor of 5-10 compared to the pre-fiare value. During 
the rapid super-fiare around MJD 55167.7 the 7-ray fiux doubles within 1 day with the 
optical and average X-ray increase of 50% and 30%, respectively. 

We find an overall correlation at all wavelengths for both long and short timescales. 
However, the unusual gamma-ray fiaring and super-fiaring activity from 4uring the period 
November-December, 2009 is not accompanied by strong emission of similar intensity in the 
optical or even in the soft X-ray bands. The hard x-ray fi ux is compa r able t o the level of 



spring 2005, e.g., ~ 3.5 x 10"^° erg crn^ s^^ in 20-200 keV 



PianetaL 



mm . This fiaring 



behavior appears to be quite different frorn othe r episodes detected in 2007 and 2008 (e.g. 



Vercellone et al. 



2009 



Donnarumma et al. 



20091 ). The synchrotron emission appears to be 



quite broad and centered around z/ ~ 10^'^ Hz. 

Striani et al. (2010) provided a first report of the AGILE-GRID data. A single 
power-law approximation gives a photon spectral index 1.66 ± 0.32 in the energy band 
0.1-1 GeV, integrating the data for two days centered at MJD 55167.7. We confirm this 
result in the analysis presented here integrating gamma-ray data for 1 day, as shown in 
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Table 1: Models parameters used for producing the spectral analysis of Fig. 2. The 
columns give the emission interval specification, the leptonic component that dominates 
the gamma-ray emission, the magnetic field B, the comoving emission radius R, the par- 
ticle energy distribution normalization K, the break energy (Lorentz factor) 7^, the min- 
imum particle energy 7mm, and the low- {(i) and high-energy spectral indices {(2) de- 
fined as A^(7) = ^ %^/[{l/lb)''^ + (7/76)^^]- The relativistic MaxweUian is defined as 
-^(7) = -f^ (7/76) exp( -7/76). 



Interval 


Model 


Component 


B 

(G) 


R 

(cm) 


K 

(cm~^) 


7b 


7min 


Ci 


C2 


Comments 


1 (pre-flare) 


two-comp. 


component- 1* 


0.6 


7 X lO^"^ 


2.2 


800 


35 


2.35 


4.5 


broken PL 


component-2 


















2 (super-flare) 


two-comp. 


component- 1 


0.6 


7 X 10I6 


2.2 


800 


35 


2.35 


4.5 


broken PL 


component-2 


0.9 


3 X lOi® 


180 


180 


1 






relativistic MaxweUian 


one-zone 


0.5 


7 X 10^^ 


2.5 


1000 


35 


2.35 


4.5 


broken PL 


3 (secondary-flare) 


two-comp. 


component- 1 


0.6 


7 X lO^'S 


2.5 


800 


45 


2.25 


4.5 


broken PL 


component-2 


0.9 


3 X lO^'^ 


170 


170 


1 






relativistic MaxweUian 


4 (post-flare) 


two-comp. 


component- 1* 


0.6 


7 X IQi® 


2.5 


800 


45 


2.25 


4.5 


broken PL 


c-()iiii)oiituil-2 



















(*)For the pre- and post-fiare intervals this set of parameters describes also the simple one- 
zone model. 
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our Fig. 2 (left panel). In addition, we found a similar spectral shape of the gamma-ray 
emission during the secondary maximum of interval-3. Our Fig. 2 (right panel) shows the 
interval-3 SED as compared with the post-fiare SED of interval-4. 

The simultaneous observations of Ijy AGILE, INTEGRAL, RossiXTE and Swift 
strongly constrain the emission models in the high-energy range. In particular, our 
interval-3 spectrum (blue solid squares of Fig. 2, right panel) represents one of the best 
constrained multifrequency spectra ever obtained for a flaring blazar from X-ray up to GeV 
energies. 

We present in Fig. 2 the results of our spectral modelling based on Synchrotron-Self 
Compton (SSC), plus contribution from external seed photons (EC). We used parameters 
similar to those alread y implemented to model pre vious gamma-ray flares of 3C 454.3. 



Vercellone et al 



fl2009[ ): 



Donnarumma et al 



(|2009[ ). We find that the pre- and post-flare 
spectra (interval-1 and interval-4) are adequately represented by a simple one-zone SSC 
model plus EC in which the accretion disk and the broad-line region provide the necessary 
soft radiation field for the inverse Compton components that dominate the X-ray through 
the GeV energies. The simple one-zone model parameters for the pre- and post-flare 
intervals are reported in Table [1] (they are also labelled component- 1 in the table). 

We can fit the super-flare with a simple one-zone model as shown in Figure |2] 
(dot-dashed curve), with the parameters reported in Table [1] Such a model requires 
(compared to the pre-flare fit parameters) an increase of the electron energy and density. 



and a slight reduction of the co moving magnetic fie 



of the blob (as also reported in 



Bonnoli et al.l 



d for the whole electron population 



20101 ). We note, however, that we used 



the optical data obtained during the rising edge of optical emission for the SED of the 
super-flare. Using only the V peak emission for the optical portion of the spectra, the 
one- zone model gives a magnetic field of 0.55 G, closer to the value for the pre-flare. 
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The time evolution of the source (e.g., the lack of strong optical emission and the 
gamma-ray spectrum during the super-flare) can be also described adopting a different 
approach. We assume a long-term rise and fall of the mass accretion rate onto the central 
black hole. This enhanced accretion causes an overall increase of the synchrotron emission 
and of the soft photon background scattered off by the primary component of accelerated 
electrons (component-1). An additional population of accelerated leptons (component-2, 
co-existing with component-1) can be introduced for the super-flare and secondary flare 
episodes. This component is a consequence of additional particle acceleration and/or 
plasmoid ejection near the jet basis. Table 1 reports the parameters that we used to 
model the Fig. 2 spectra (the two-component models are reported as thick solid lines for 
the super-flare and for the secondary flare). For all intervals, we assume the presence 
of component-1. This is the only component acting for the pre- and post-flare intervals, 
with the parameters given in Table 1. The super-flare (secondary-flare) parameters for the 
component-1 are the same as the pre-flare (post-flare) interval. We assumed a bulk Lorentz 
factor r = 25, a jet angle with respect with the hue of sight 9 — 1.2°, an accretion disk of 
bolometric luminosity = 6 x 10^^ergs~^slowly decaying toward = 5 x lO^^ergs^^. A 
broad line region located 0.5 pc from the black hole reflects 5% of the disk power toward 
the emitting regions. 

The component-2 energy distribution that better reproduces our gamma-ray spectral 
data is a relativistic MaxweUian of characteristic energy 7^ ~ 180. Interestingly, this 
component appears to be strongly energized but not yet modifled by additional non-thermal 
acceleration. 

This two-component approach avoids the problem of explaining the time variability 
of the physical parameters of the whole electron population of the one- zone model, at the 
price of adding an additional varying component superimposed with the pre-existing one. 
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To summarize, our multifrequency data for the December, 2009 flare of provide a 
wealth of very important information on this puzzling and fascinating blazar. We find 
that is characterized by strong broad-band spectral variability, and that the modelling of 
the peak gamma-ray emission episodes suggests more elaborate models than the standard 
one-zone SSC-I-EC models of bright blazars. 
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